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Abstract. The soon-to-be-launched Space Infrared Telescope Facility 
(SIRTF) shall produce image data with an a-posteriori pointing knowl¬ 
edge of 1.4" (lcr radial) with a goal of 1.2" in the International Celestial 
Reference System (ICRS). To perform robust image coaddition, mosaic 
generation, extraction and position determination of faint sources, the 
pointing will need to be refined to better than a few-tenths of an arc- 
second. This paper summarizes the pointing-refinement algorithm and 
presents the results of testing on simulated data. 


1. Introduction 

One of the goals of astronomical image data acquisition is to infer the pointing in 
an absolute coordinate system as accurately as possible, in this case the celestial 
reference system. Instabilities in telescope pointing and tracking however inhibit 
us from achieving this goal. SIRTF for instance, is expected to provide pointing 
and control of at least 5" absolute accuracy with 0.3" stability over 200 sec (lc 
radial). In the ICRS, the star-tracker assembly provides an a-posteriori pointing 
knowledge of 1.4". The end-to-end pointing accuracy is a function of the inherent 
star-tracker accuracy, the spacecraft control system, how well the star-tracker 
bore-sight is known in the focal plane array (instrument) frame, and variations 
in the latter due to thermo-mechanical deflections. The SIRTF observatory will 
have a data flow rate of > 10, 000 science images per day. This will require an 
automated, self-consistent means of refining the celestial pointing as robustly as 
possible. 

The conventional method to refine the pointing is to make comparisons of 
astrometric sources with positions known to better than a few percent of the 
observed positional errors. The primary motivation for pointing refinement is 
to enable robust coaddition of image frames in a common reference frame so 
that source extraction and position determination to faint flux levels can be 
performed. Moreover, refinement in an absolute (celestial) reference frame will 
enable robust cross-identification of extracted sources with other catalogs. 

For these purposes, we have developed a stand-alone software package u point- 
ingrefine ” with a goal to generate science products with sub-arcsecond pointing 
accuracy in the ICRS. The software can refine the pointings and orientations of 
SIRTF images in either a “relative” sense where pointings become fixed relative 
to a single image of a mosaic, or, in an “absolute” sense (in the ICRS) when 
absolute point source information is known. 
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Figure 1. A simple three image mosaic. 


2. Software Summary 

As part of routine pipeline operations at the SIRTF Science Center (SSC), all 
input images are pre-processed for instrument artifact removal and pointing 
data attached to raw FITS images. Following this, point source extraction 
is performed on individual frames for input into pointingrefine. The software 
expects point source lists adhering to the format produced by the SSC source 
extractor. The pointingrefine software performs the following: 

1. Reads source extraction tables, FITS images from input lists and if abso¬ 
lute refinement is desired, a list of absolute point source positions. 

2. Point source position and flux matching is performed between all possi¬ 
ble image pairs in the input list. This includes absolute point sources if 
available. 

3. Transform correlated point source positions and uncertainties to a Carte¬ 
sian fiducial (mosaic) reference frame. 

4. Set up global minimization equations involving relative offsets between all 
correlated point source positions. 

5. Solve a linear matrix equation for translational and rotational offsets in 
mosaic reference frame for each input image and compute the full error- 
covariance matrix. 

6. Apply these offsets to the pointing centers (x c , y c ) of all images. Transform 
to RA, Dec to obtain refined celestial pointings and orientations. 

7. Results output: Write new refined pointing keywords to FITS headers and 
to a table file with diagnostic information. 


3. Algorithm 

A brief outline of the refinement algorithm (primarily steps 4, 5, and 6 above) is 
as follows. Consider the simple three image mosaic in Figure 1. Image “1” defines 
the “fiducial” reference frame. The circles represent point sources detected from 
each overlapping image pair transformed into the fiducial frame. The filled circles 
are sources extracted from image n and the open circles are sources extracted 
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from either image 1 or m. The correlated source pairs are slightly offset from 
each other to mimic the presence of pointing uncertainty in each raw input 
image. The Cartesian coordinates of a correlated point source common to an 
image pair are related by: 


= x?-{y™-y™)50 m + 5X m 
y? -> y? = y r-(xT-x™)66 m + 8Y m , 


( 1 ) 


where 59 m , 5X m , 5Y m are rotational and Cartesian offsets respectively in the 
fiducial mosaic frame. We have made the approximation sin 59 ~ 59 (cos 59 ~ 1) 
since uncertainties in measured position angles are expected to be small (< 20 // ). 

We define a cost function L, representing the sum of the squares of the 
“corrected” differences of all correlated point source positions in all overlapping 
image pairs (m,n): 


where 




( 2 ) 


Aaf’ n = a 2 (x?) + a\x?) (3) 

a yT n = * 2 (yD + °\y?) 


and the a 2 represent variances in extracted point source positions. The function 
L (Equation 2) is minimized with respect to the Cartesian offsets ( 59 m , 5X m , 
5Y m ) for an image m which contains point sources in common with all other 
images n. At the global minimum of L, the following conditions hold: 


8L _n. dL _n. dL 

859 m ~ °’ 85X m ~ °’ 85Y m 


( 4 ) 


Evaluating these partial derivatives leads to a set of three simultaneous equations 
for each image in our mosaic. In general, for N correlated images, we will have 
3 (IV — 1) simultaneous equations in 3(IV — 1) unknowns. It is U N — 1” because 
we exclude the reference image frame which by definition has the constraint 
89 m = 5X m = 5Y m = 0. The 3 (N — 1) system of equations can be solved 
using linear sparse matrix methods since a large number of the matrix elements 
will be zero. We use the UMFPACK 1 library which is adapted for solving large 
unsymmetric matrix systems. 

When absolute astrometric references are available, the “fiducial” reference 
image (mosaic) frame is treated like a single input image, which contains the 
absolute point source positions. When the input images become refined relative 
to this fiducial image, they in reality become absolutely refined (in the ICRS). 
The presence of absolute point sources also reduces the effect of “random walks” 


1 http:// wot. cise.uf1.edu/research/sparse/umfpack/ 
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Figure 2. Distribution of radial offsets between “true” (noiseless) and 
simulated pointings (with noise) - broad histogram, and between true 
and refined pointings - narrow histogram. 


in offset uncertainties with distance if a single input image were chosen as the 
reference instead. Once Cartesian offsets (58 m , 5X m , SY m ) in the reference 
(mosaic) frame are computed, the pointing centers (x c , y c ) are corrected (by 
use of Equation 1) and transformed back to the sky to yield refined pointings. 
Image orientations are refined in a similar manner, but in this case, we need 
to transform at least two fiducial points per image to uniquely determine the 
orientation. 


4. Simulations 

The Infrared Array Camera (IRAC) on SIRTF will perform simultaneous imag¬ 
ing at four bands spanning the range « 3.6 fjam to 8/rni. Each array consists 
of 256 x 256 1.2" pixels. We simulated a mosaic of 800 IRAC (3.6 pm) “truth” 
images (i.e. with no pointing error) with each image containing randomly dis¬ 
tributed point sources. Input image overlap coverage was ~ 50%. A second 
set of 800 images was simulated with random errors added to the pointing key¬ 
words of image headers. The errors were drawn from a Gaussian distribution 
with mean radial error (S') ~ 1.4". An absolute point source list was also sim¬ 
ulated by extracting the brightest sources with smallest centroiding errors from 
each “truth” image. The average number of “absolutes” per input image was 
10. The SSC point source extractor was then used to extract ~ 50 point sources 
from each input image (with pointing error). 

Figure 2 shows the results of our simulation where we compare the distribu¬ 
tions of radial offsets relative to “truth” pointings before and after refinement. 
The refinement is better than 85% for almost every image. The main limitation 
is full knowledge of the Point Response Function (PRF) to reduce centroiding er¬ 
rors in source extractions. However, we expect extraction centroids better than 
0.1" with better sampled PRFs. This will give us the sub-arcsecond absolute 
pointing accuracy sought in SIRTF’s imaging detectors. 
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